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THEORETICAL  COMPUTATIONS  OF  SOUND  REFLECTION 
FROM  A LAYERED  OCEAN  BOTTOM 


INTRODUCTION 

The  U.  S.  Naval  Oceanographic  Office  is  engaged  in  o program  to  de- 
termine the  influence  of  bottom  composition,  structure,  and  topography  on 
sound  reflection  from  the  ocean  bottom.  As  one  step  towards  accomplishing 
this  objective,  a computer  program  was  prepared  to  enable  theoretical  investi- 
gations of  the  effect  of  layered  sediments  on  sound  reflection-  The  computer 
mathematical  mode!  permits  detailed  and  systematic  calculations  of  bottom  re- 
flection coefficients  and  bottom  reflection  loss  as  the  following  factors  are  varied 
parametrically:  number  of  reflecting  layers,  layer  thickness,  sediment  sound 
velocity,  density,  sediment  absorption,  and  frequency.  This  paper  will  illustrate 
how  a layered  ocean  bottom  model  can  be  used  to  develop  an  understanding  of 
ocean  bottom  reflectivity. 

DESCRIPTION  OF  COMPUTER  PROGRAM 

The  computer  program  computes  the  complex  reflection  coefficient,  phase 
shift,  and  reflection  loss  for  a plane  wave  obliquely  incident  on  a layered  ab- 
sorbing bottom.  Expressions  for  the  reflection  coefficient  are  derived  from  the 
wave  equation  by  utilizing  steady  state  solutions  to  the  differential  equation  for  a 
"Voigt  Liquid."  Matrix  algebraic  operations  after  the  method  of  Haskell  (1)  are 
employed . 

The  mathematical  derivation  assumes  a multilayered  ocean  bottom  in  which 
the  sediment  layer  thickness  (D),  sound  velocity  (C),  density  ( P),  and  absorption 
(a)  are  assigned  arbitrary  values.  The  number  of  layers  considered  is  not  limited 
and  the  lowermost  layer  in  the  sedimentary  sequence  is  considered  to  be  semi- 
infinite. The  sedimentary  layers  are  assumed  to  be  fluid;  that  Is,  the  shear  mod- 
ulus |X  - 0.  The  layers  are  further  assumed  to  be  absorbing  and  to  have  plane 
parallel  interfaces.  The  pressure  and  particle  velocity  one  continuous  across  the 
interfaces  between  the  layers. 

The  computer  program  was  developed  and  written  by  Dr.  James  Dorman,  as 
a consultant  to  Alpine  Geophysical  Associates  (2),  under  contract  to  the  U»  S, 

Naval  Oceanographic  Office.  A detailed  discussion  of  the  theory  involved  in  ] , 

the  development  of  the  program  can  be  found  in  reference  (2) . The  following  » . ' 

equations  are  presented  as  background  information  as  the  fundamental  relation-  j-" n 

ships  utilized  In  developing  the  layered  model. 
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The  velocity  potential  of  compressicnol  wove*  that  travel  In  Hie  fluid 
must  satisfy  Hie  differential  equation  for  a "Voigt  Liquid, " 
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A solution  of  the  differential  equation  fs 
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where  the  terms  of  d>  containing  A and  A'  represent  oblique  downgoing  and 
upgolng  waves,  respectively,  and  x ± z Kxn  e represents  the  horizontal  dls— 
tosice  the  wave  front  travels  in  time  t. 


Substituting  equation  2 in  equation  lb  leads  to 
u?z  $ ( a*+ 1 w a'* ) ( i + tan3® ). 


(3) 


This  equation  is  satisfied  by  complex  f . Substituting  F~Hr  * Z 
3 yields 


In  equation 
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where  the  negative  root  of  ' it  taken  to  represent  attenuation  in  the  direction 
of  propagation  and  the  positive  root  of  ij  represents  the  horizontal  wave  number . 
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Tite  expression  for  the  velocity  potential  it  used  to  write  a matrix  to  solve 
for  the  downward  particle  velocity  and  fluid  pressure.  By  matrix  inversion  an 
expression  for  a fluid  layer  with  Voigt-type  attenuation  Is  arrived  at,  and  by 
means  of  recurrent  matrices  the  reflection  coefficient  for  an  acoustic  wave  re- 
flected from  n layers  is  determined. 

The  computer  program  was  written  in  Fortran  II  for  an  IBM  7074  computer. 

There  are  two  types  of  input  cards  for  the  program.  The  first  type  specifies  the 
four  parameters  of  each  layer  and  one  card  per  layer  is  used.  The  parameters 
specified  are  layer  thickness,  layer  sound  velocity,  density,  and  absorption. 

The  absorption  constant  is  specified  as  on  amplitude  per  unit  length,  such  as 
decibels  per  meter. 

The  second  type  of  card  is  a frequency  card  which  specifies  the  frequency  ?n 
cycles  per  second. 

All  layer  cards  occurring  in  a continuous  sequence  are  considered  to  be  layers 
of  a single  model . The  first  card  is  treated  as  the  first  layer  below  the  ocean  bot- 
tom and  the  last  card  is  treated  as  a semi-infinite  half  space.  The  sound  velocity 
and  density  of  the  water  layer  are  fixed  in  the  program  but  may  be  readily  changed 
to  any  values. 

The  output  consists  of  completely  annotated  tables  of  reflection  coefficients, 
phase  shifts,  and  reflection  loss  at  grazing  angle  intervals  of  1 degree.  This  in- 
terval may  be  easily  changed  as  desired. 

TYPICAL  COMPUTED  CURVES 

All  of  the  computed  curves  illustrated  in  Figures  1 to  II  were  arrived  at  by 
assuming  or  inventing  typical  ocean  bottom  conditions.  The  values  used  for  sedi- 
ment sound  velocity  are  consistent  with  the  data  of  Hamilton  et  ol.  (3),  Shorn  way 
(4),  and  Sutton  et  al , (5)  , This  was  done  to  facilitate  the  computations  and  to 
permit  each  of  the  specified  variables  to  be  varied  parametrically.  By  so  doing 
a systematic  investigation  of  the  effect  of  frequency,  sediment  sound  velocity, 
layer  thickness,  absorption,  and  the  number  of  reflecting  layers  on  sound  reflection 
from  the  ocean  bottom  can  be  performed.  While  it  is  apparent  that  many  com- 
binations of  there  variables  may  exist  it  is  possible,  at  this  time,  to  present  only 
a limited  number.  Each  of  the  figures  presents  curves  of  bottom  reflection  loss  as 
a function  of  grazing  angle  for  a given  set  of  conditions.  The  following  units 
should  be  used  when  referring  to  the  illustrations:  sound  velocity  (C)  in  feet  per 
second,  density  (f3)  in  grams  per  cubic  centimeter,  layer  thickness  (D)  in  feet, 
and  absorption  (a)  in  decibels  per  foot. 
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Figure  1 to  6 represent  computed  I mi  values  for  a three  layer  model  in 
which  the  thickness  of  the  second  layer  w os  arbitrarily  specified  as  5 feet.  In 
addition/  for  this  model  the  bottom  was  assumed  to  be  nonabsorbing  (a  = 0).  in 
all  coses  the  first  layer  specified  (layer  1)  represents  the  semi-infinite  water 
layer  overlying  the  first  sediment  layer  (foyer  2), 

The  dependence  of  bottom  reflection  loss  on  frequency  for  a layered  ocean 
bottom  can  be  $t«n  in  Figures  1 and  2.  A comparison  of  the  i,  4,  and  12  Icc 
traces  in  these  figures  shows  thct  os  frequency  Increases  reflection  loss  becomes 
more  dependent  on  grazing  anpje.  The  observed  dependence  results  from  inter- 
ference effects  produced  by  reflections  from  the  water-sediment  interface  ami 
from  the  subbottom  layer,  and  is  seen  to  occur  as  fluctuations  In  the  reflection 
loss  curve  os  a function  of  grazing  angle.  In  Figure  1 the  1 kc  curve  appears  t s 
a smoothly  oscillating  curve  with  2 peaks;  the  4 Icc  curve  has  6 peaks;  and  the  12 
kc  curve  Is  rapidly  oscillating  with  18  peaks.  The  variations  in  reflection  loss 
with  grazing  angle  range  from  a low  of  5 decibels  at  1 kc  to  a high  of  about  18 
decibels  at  12  kc.  Since  the  bottom  was  specified  to  be  nonabsorbing,  total  re- 
flection occurs  at  the  critical  angle.  Figure  1 represents  a model  where  Is 
the  maximum  velocity  and  the  minimum  velocity  (C^XL^C-)'  w^e  ‘n 
Figure  2 CL  U the  maximum  velocity  and  C_  the  minimum  The  some 

frequency  dependence  illustrated  In  Figure!  is  seen  in  Figure  2,  but  the  amplitude 
of  the  fluctuations  Is  nor  as  forge.  The  maximum  variation  is  about  8 decibels 
and  at  a grazing  angle  of  about  35  degrees  the  oscillations  are  almost  completely 
dampened.  These  figures  indicate  that  o layered  ocean  bottom  can  explain  fluctua- 
tions In  bottom  reflection  lass  as  a function  of  grazing  angle  and  that,  for  a non- 
absorbing bottom,  reflection  loss  data  need  not  exhibit  a uniform  dependence  on 
frequency. 

The  effect  of  sound  velocity  variations  in  layer  3 (CL)  can  be  estimated  from 
an  examination  of  Figures  3 and  4.  Except  for  shifts  In  tne  critical  angle  the  shape 
of  the  curves  remains  relatively  the  same  as  Increases;  however,  there  are  notice- 
able differences  in  the  amplitude  of  the  fluctuations  that  appear  to  be  controlled  by 
CL.  it  can  be  seen  in  Figure  3 that  for  angles  less  than  30  degrees  the  reflection 
Idas  Increases  for  decreasing  values  of  C^.  A transition  zone  is  observed  at  30  de- 
grees. For  anglos  greater  than  this  the  amplitude  of  the  oscillations  increase  as 
C Increases  and  the  maximum  and  minimum  loss  values,  in  this  region,  are  associated 
with  Increasing  CL . The  smooth  curve  in  this  figure  represents  o two  foyer  model 
formed  by  foyers  fond  3 for  CL  = 5200  feet/ second.  It  Is  apparent  that  the  two  k/er 
model  determines  the  minimumioss  or  lower  limit  of  the  three  layer  case.  If  should 
be  noted  that  Cj  is  the  minimum  velocity  in  Figuro  3,  whereas  C.  is  the  minimum  In 
Figure  4.  Figure  4 mere  dramatically  illustrate*  the  increase  in  tne  amplitude  of  the 
fluctuations  as  increases. 
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FIGURE  4.  EFFECT  OF  C ON  BOTTOM  LOSS  FOR  C,  MINIMUM 
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figures  5 and  6 demonstrate  the  importance  of  layer  2 sound  velocity  ( C ^ 
in  computing  reflection  loss.  In  Figure  5 C^Cy  but  P^Z,  <P  C . A,C2.n- 
creases  the  number  of  oscillations  increase,  while  the  amplitude0 decreases.  In 
the  next  illustration  (Figure  6)  the  maximum  velocity  is  determined  by  C2  and 
the  same  dampening  effect  observed  in  Figure  2 can  be  seen  again.  This  smoothing 
effect  is  soen  to  shift  towards  the  higher  grazing  angles  as  C2  increases.  As  in 
Figure  5 the  amplitude  of  the  fluctuations  decreases  as  C2  increases 

The  influence  of  layer  2 tl  'ckness,  in  a three  layer  model,  is  shown  in  Figure 
7.  The  layer  thickness  was  varied  with  respect  to  wavelength.  As  the  thickness 
increases  from  a quarter  wavelength  to  8 wavelengths  the  reflectivity  curve  changes 
from  a smooth  curve  to  one  having  10  peaks.  As  expected,  for  notmal  incidence 
(90  degrees),  a maximum  loss  value  is  obtained  for  the  quarter  wavelength  and  o 
minimum  value  for  the  half  wavelength  layer. 

Figure  8 illustrates  the  result  of  increasing  the  number  of  reflecting  layers  to 
five.  As  was  fho  case  in  Figures  1 and  2 the  fluctuations  increase  os  the  frequency 
increases;  however,  in  this  case  the  reflectivity  curves  are  not  nearly  as  orderly. 

This  results  from  the  interference  effects  produced  by  four  reflecting  Interfaces  as 
the  grazing  angle  varies. 

To  illustrate  the  effect  of  absorption  in  the  sediments  on  the  bottom  reflection 
curves,  Shumway's  (4)  20  to  40  kc  absorption  data  were  extrapolated  according  to 
a square  root,  first  power,  and  square  dependence  on  frequency . The  results  of 
this  test  are  illustrated  in  Figures  9,  10,  ond  11.  For  the  case  of  reflection  with 
absorption  there  is  no  critical  angle  and  the  reflection  never  becomes  total.  This 
is  illustrated  in  Figure  9 for  1 kc.  Absorption  not  only  eliminates  the  critical  angle 
and  Increases  the  loss  values  for  angles  less  than  30  degrees,  but  also  tends  tc  dampen 
the  fluctuations.  This  is  most  noticeable  for  the  case  of  a 0.261  decibel  per  foot, 
where  the  peak  loss  at  45  degrees  of  about  26  decibels  for  a - 0 deceases  to  about 
19  decibels.  As  the  absorption  decreases  this  effect  becomes  less  noticeable.  Figure 
10  depicts  this  same  effect  at  4 kc.  Notice  the  extreme  dampening  for  a ~ 1 ,01 
decibels  per  foot.  For  a “ 0. 108  an  additional  peak  has  formed  in  the  region  less 
than  the  critical  angle.  Figure  11  represents  a situation  where  D was  Increased  to 
about  20  feet  as  compared  to  5 feet  in  Figure  10.  The  number  of  oscillations  has 
increased  and  the  effect  of  absorption  is  noticeable  by  the  dampened  peaks.  It  ap- 
pears that  absorption  tends  to  eliminate  some  of  the  effects  of  constructive  and  de- 
structive Interference,  particularly  far  larger  values  of  absorption. 
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FIGURE  5.  DEPENDENCE,  OF  BOTTOM  LOSS  ON  C„. 
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FIGURE  6.  DEPENDENCE  OF  BOTTOM  LOSS  ON  C 
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BOTTOM  LOSS  AS  A FUNCTION  OF  GRAZING  ANGLE  FOR  VARYING  LAYER 
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GRAZING  ANGLE  (DEGREES) 

FIGURE  9.  EFFECT  OF  ABSORPTION  AT  1 KC  ON  BOTTOM  LOSS. 
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COMPARISON  WITH  EXPERIMENTAL  DATA 

in  on  attempt  to  establish  the  usefulness  of  the  layered  model  for  predicting 
reflection  Iocs  a comparison  was  made  of  measured  loss  values  obtained  in  a water 
depth  of  about  2800  fathoms,  in  the  Hatteras  Abyssal  Plain,  with  values  computed 
from  a short  sediment  core.  A 2.5  foot  core  was  taken  at  the  location  of  a soft  of 
6 kc  bottom  loss  measurements  made  with  explosive  sound  sources.  The  sediment 
grain  size,  porosity,  density  and  various  other  physical  properties  ‘..  ere  measured 
directly  from  the  core.  The  Shumway  semi-empirical  relationship  between  porosity 
and  velocity  was  used  to  compute  sediment  sound  velocities.  As  indicated  by 
Hamilton  (6)  these  velocities  were  then  corrected  for  the  differences  in  temperature 
and  pressure  between  S hum  way's  laboratory  measurements  and  in  situ  conditions. 
Wilson's  (7)  tables  for  the  speed  of  sound  in  sea  water  were  used  to  moke  this 
correction.  Absorption  values  were  extrapolated  from  Shumwoy's  high  frequency 
data  according  to  the  first  power  frequency  relation  used  by  Mackenzie  (8)  and 
Ccle  (9).  The  measured  density  and  estimated  acoustic  properties  were  then  used 
to  establish  five  thin  layers  for  use  in  the  ir  'el . 

A comparison  of  the  predicted  awa  ;.:ured  loss  values  at  6 kc  is  presented  in 
Figure  12.  Except  at  a grazing  angle  of  70  degrees  the  overall  agreement  appears 
to  be  good.  The  mea«»red  values  represent  the  reflection  loss  at  the  water-sediment 
interface  and  the  short  core  used  for  the  computations  should  be  adequate  due  to 
the  extremely  short  pulse  of  the  explosive  sound  source. 

A similar  comparison  of  computed  and  measured  losses  was  made  for  a set  of 
3.5  kc  measurements  made  in  about  2800  fathoms  of  water  in  the  Sohm  Abyssal 
Plain.  A short  cere  was  also  used  for  the  computations  and  the  results  are  presented 
in  Figure  13.  The  overall  agreement  between  the  computed  and  measured  loss 
values  appears  to  be  quite  good. 

Other  unpublished  comparisons  do  not  always  exhibit  such  agreement;  how- 
ever, this  does  not  necessarily  suggest  a deficiency  in  the  theory  or  model,  but 
may  very  well  indicate  that  not  enough  is  known  about  the  exact  relationships 
between  the  acoustic  and  physical  properties  of  bottom  sediments.  While  the 
studies  presented  in  references  4,  5,  and  6 are  excellent  studies  one  has  only  to 
compare  the  results  of  each  to  see  that  differences  exist  between  the  various  en- 
vironments sampled  and  that  a universal  relationship  between  the  acoustic  and 
physical  properties  of  bottom  sediments  probably  does  not  exist. 
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FIGURE  12.  COMPARISON  OF  COMPUTED  AND  MEASURED 
REFLECTION  LOSS  AT  6 KC. 
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FIGURE  13.  COMPARISON  OF  COMPUTED  AND  MEASURED 
REFLECTION  LOSS  AT  3.5  KC. 
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A computer  mod©!  has  been  developed  to  permit  systematic  theoretical  in- 
vestigations to  aid  in  the  interpretation  and  prediction  of  ocean  bottom  reflec- 
tivity, In  addition,-  such  a mode!  can  be  of  value  in  planning  acoustic  measurement 
progress.  The  layered  ocean  bottom  mode!  can  be  used  to  explain  the  effect  of 
frequency,  layering,  and  sediment  properties  on  bottom  reflection  loss  as  a function 
of  grazing  angle. 


Theoretical  computations  cGn  be  used  for  interpreting  acoustic  data  for  the 
purpose  of  defining  the  acoustic  properties  of  bottom  sediments;  however,  extreme 
cere  must  be  exercised  v/hen  so  doing,  as  It  is  apparent  that  various  combinations 
of  sound  velocity,  density,  absorption,  and  layer  thickness  can  produce  the  same 
reflectivity  curve. 

Limited  comparisons  of  computed  and  measured  reflectivity  In  abyssal  plain 
regions  have  shown  that  experimental  data  ore  consistent  with  a model  that  con- 
siders the  ocean  bottom  to  ba  flat,  absorbing,  layered,  and  fluid.  The  usefulness 
of  the  model  for  predicting  bottom  reflectivity  is  limited  by  our  knowledge  of  the 
exact  properties  of  the  ocean  bottom.  ^ 
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